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By O. Zeman ' 

1. Motivations and objectives 

The purpose of this work was to investigate the turbulent constitutive relation 
when turbulence is subjected to solid body rotation. Rotating turbulent flows exist 
in many industrial and geo- and astrophysical applications. 

2. Accomplishments 

2.1 Note on spectra and decay of rotating homogeneous turbulence 

Recently, Squires, Chasnov, Mansour, & Cambon (1993) (hereon SCMC) ad- 
dressed the problem of asymptotic behavior of homogeneous turbulence. Briefly, to 
summarize their results and findings, they applied a spectral LES method to achieve 
an asymptotic, equilibrium evolution of initially isotropic turbulence subjected to 
solid body rotation of angular speed ft. The computations, which were run for 
times of order O(10 3 ) of initial turbulence time scales, confirmed the prediction of 
the asymptotic decay laws. The authors used two different (initial) spectral forms 
of the energy spectrum E(k) at low wavenumber k: 

E(k) = 2irk 2 A 0 + .... and E(k) = 2nk A A 2 + .... 

In nonrotating turbulence the two spectral forms are known to produce different 
time decay exponents (n) of the turbulent kinetic energy ^q 2 = / 0 E\k)dk oc t 
For the k 2 spectrum, n = 6/5, and for the k* spectrum, n = 10/7. In the presence 
of rotation the following asymptotic decay laws were proposed in SCMC: 

q 2 oc A 2/5 t -3/5 ft 3 ^ 5 ( k 2 spectrum) (1) 

q 2 oc Al /7 r 5 ' 7 n 5 ' 7 ( k 4 spectrum) (2) 

The above laws were confirmed by the LES computations within a few percent. 
Computations also indicated that the rotating turbulence has a tendency toward 
a two-dimensional state in the sense that the spectral energy tends to concentrate 
at wavenumbers normal to the rotation axis, i.e. the gradients with respect to the 
wavenumber parallel to ft become relatively small. At the same time the turbulence 
remained remarkably close to isotropy if measured by departure from the isotropy 
tensor bij =< u,Uj > jq 2 - 1/36, r This suggests a turbulence structure consisting 
of vortices aligned with the rotation axis and of jet-like (fluctuating) flow parallel 
to the rotation axis. 

The purpose of this note is to explain the behavior of the rotating turbulence on 
the basis of a model for the spectral energy transfer, and to propose modification 
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of the turbulence spectrum when the rotation is much more rapid than the eddy 
turnover time scale. We start with a simplified Lagrangian description of the relation 
between stress and mean strain in rotating turbulence. The resulting relationship is 
then used to describe the energy transfer from larger to smaller scales in the spirit 
of the analysis described in Tennekes & Lumley (1972) for nonrotating turbulence. 

Denoting Lagrangian fluctuating velocity components by Vi and taking the ro- 
tation vector to be O = (0,0, fl), one can write equations for the motion of fluid 
particles originating at some point in space and time (ao,f Q ) as follows: 

dvi 

— = -VjUij + 2 tijsQvj + n i (3) 


Here, J7, represents a background mean velocity field which is considered as slow 
varying with respect to the characteristic Lagrangian (turbulence) time scale 77, 
and length scale (oc 977,). H\s are random forcing terms comprising the effect of 
the fluctuating pressure and higher order correlations. The viscous terms are taken 
to be negligible on account of the high turbulence Reynolds number assumption 
(Rer oc q 2 ri/u >> 1). The velocities 17 are functions of position and time v,(X,t) 
of the fluid particle, with the initial position X(f G ) = ao. Because turbulence 
is statistically homogeneous, we shall suppress the space dependence and utilize 
the ensemble-average identity < ViVj > (t) =< U{Uj > (t), i.e., the one-point 
Lagrangian averages (over all initial locations) are equal to the Eulerian averages 
(over the flow volume). A useful reference for Lagrangian description of turbulence 
is Monin &: Yaglom (1971). 

Now, neglecting the effect of IT, and assuming that the gradients {7 l 3 of the 
slow- varying velocity field (in the rotation direction) are negligible, it is possible to 
formulate a stress-strain relation < V\V 2 >oc — S12 ( Sij = + Uj,i ) is the slow 

strain tensor). This is achieved first by integrating (3) to obtain expressions for 
v u v 2 : ^ ^ 

v\(t) = v\o(t 0 ) - Uij I Vj(t , )dt t + 2fl f v 2 (t')dt' 

Jt 0 Jt 0 

v 2 (t) = v 2o (t 0 ) - U 2 ,j f t )j(t')dt 9 — 2ft f vi(t f )dt' 

Jto Jto 

Further manipulations yield an expression 


< v x v 2 >—< u\u 2 >= — [a 


W r i 


1 + Ac 2 tIQ(2Q — Ri 2 ) 


]2Si2 = —2vj'S\ 2 . 


(4) 


Here, all the required Lagrangian time scales such as Tl q =< >~ l J 0 °° < 

i; Q (^)t; a (f 0 ) > dt* were written for simplicity as a single time scale 77,, which is in 
turn proportional to the turbulence time scale (,/u 1 (£ and u f being chairacteristic 
length and velocity scales); f?,j is the asymmetric complement to Sij. The presumed 
differences in timescale magnitudes axe absorbed in the free constants C { . Clearly, (4) 
expresses a turbulent constitutive relation in the presence of rotation; the effective 
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eddy viscosity ut is evidently diminished by a factor depending on (17, ft) 2 . Although 
in the following analysis the flow field represented by Sij and Rij is taken as random, 
it is of interest to interpret (4) considering a homogeneous rotating shear flow with 
5j 2 = R 12 = \U \ t 2 = |S. Then we note that the nature of the constitutive relation 
(4) is such that ut is maximized for (ft/ S) ma x = 0.125. On the other hand, the 
LES results of Bardina et al. (1985) and linear stability analysis (see e.g. Speziale, 
1991) indicate the maximum turbulence amplification at ($l/S) ma z = 0.25. Since 
in the following analysis R i2 is neglected, the exact value of (ft/5) ma i is irrelevant 
to our problem. A stress-strain relation analogous to (4) has also been derived in 
cylindrical coordinates for a turbulent line vortex by Zeman (1994a) (comparable 
when the vortex flow is in solid body rotation i.e. when Uazimuth = r ft)- 

Now we shall relax the relation in (4) so that S12 = S(k') represents the strain 
of eddies (of size oc 1 /k') larger than the wavenumber k of the stress < tii«i 2 > (AO- 
Following the line of reasoning in Tennekes k Lumley (Section 8.4) concerning the 
spectral transfer in nearly isotropic turbulence, the spectral energy flux T(A) across 
the wavenumber k in the inertial subrange is effected mainly by local interaction 
so that T(A) oc -S tJ (k') < UiU^k") >oc i/ t (A")S 2 (A') where, approximately, A/3 < 
k' < k and k" = 3k'. As shown in Tennekes k Lumley, the quantities at k' or k" 
are directly related to the same quantities at A, thus e.g., S(k') oc ( E(k)k 3 ) 1 / 2 oc 
r[ l (k"). Utilizing (4) to express T(k) in terms of quantities (depending now on 
E(k), k , ft) and neglecting the contribution ft_R 12 in (4) (Rn{k) is a random 
quantity with zero mean and |.ff 12 | << ft), we obtain a relationship 


T(A) oc u T {k")S 2 (k') oc 


(Ekf/ 2 k 


In the inertial subrange the spectral energy flux T(A) across each wavenumber is 
constant and equal to the dissipation e, and the above equation can be written as 


e 


— 3 /2 (Ek)^k 

, . „ ft 2 ’ 
1 + c 3 #pr 


(5) 


where a is the Kolmogorov constant and C3 is another free coefficient. Evidently, 
(5) represents an implicit relation for the energy spectrum E(k,e,£l) = 0 in the 
presence of rotation. For more insight into the meaning of (5), it is useful to define 
a rotation (cut-off) wavenumber fcn 


k a = (^-)'l\ (6) 

which delimits the region of the spectrum where the rotation effects are impor- 
tant, i.e., k < kit; (note that A^ 1 is analogous to the Ozmidov length in stratified 
turbulence). In the region where k << An, (5) results in an explicit expression 

E(k) oc e 2 / 5 ft 4 / 5 A- n / 5 , 


(7) 
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FIGURE 1. Sketch of the turbulence energy spectrum subjected to rotation; k Q 
and k n are, respectively, the large-eddy and Kolmogorov scale wavenumbers, kQ is 
defined in (6). 


while for k >> kQ the Kolmogorov inertial subrange E(k) = ae 2 / 3 fc“ 5 / 3 is recov- 
ered. The sketch of the spectrum (with exaggerated slope change) is depicted in 
Fig. 1. A general solution of (5) (with C3 = 1) is obtained in the form 


x 


= ( 


1 + 2 y + y 2 


) 1/4 > 


( 8 ) 


where x = k/kQ and y = Ek 3 /U 2 . Fig. 2 shows the solutions of the above equation 
emphasizing the rot at ion- affected range by plotting Ek 11 / 5 and Ek 5 / 3 . It is seen 
that the spectrum of the form (7) is approximately valid for k/kQ < lO” 1 . It should 
be noted that expressing the eddy viscosity (in square brackets) in (4) in terms of 
spectral quantities at a given k, one obtains 


u T (k) = 


(. E/k ) 1 ! 2 
1 + c 4 /y{k)' 


( 9 ) 


hence the parameter fl 2 /Ek 3 is the measure of the damping effect of rotation on 
the local eddy viscosity. It is of interest that the same parameter appears in the 
expression for the subgrid-scale eddy viscosity in the LES of SCMC. Although the 
functional dependence of 1 'r(k) on y is far more complicated, both CSMC and 
expression (9) give the same asymptotic dependence ur(k) cxy if y « 1, 

If there exists a self-similar spectrum as sketched in Fig. 1 (with k 0 < kQ << 
fc^), then (7) also contains information concerning the turbulence energy decay. 
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Figure 2. Spectral energy quantities Ek 11/S ( ) and Efc 5/3 ( ) vs. 

1/jfcjj, baaed on the solution to (8). par 

Assuming a k 4 spectrum for k < k a , then oc A 2 (fc®),t oc — € » us ' n g (7) one 

° bt “ nS (I") 

The exponent n = 10/21 is lower than n = 5/7 in (2) derived from the dimensional 
analysis in SCMC. If, however, the rotation damping factor in (5) can be generalized 
to 1 + c 4 y _1 oc y~ m , one obtains the relation n = 10/(7 + 14m). Hence, to satisfy 
the decay exponent n = 5/7 proposed in (2), m must take on the value m = 1/2. 
This leads to a spectral form 

E(k) oc(eQ) 1/2 k~ 2 , 


different from (7). Analogous relations can be obtained for the k 2 spectrum to 

In summary, from Lagrangian analysis a relation between turbulent stress and 
strain in rotating homogeneous turbulence was inferred. This relation was utilized 
to derive the spectral energy flux and, ultimately, the energy spectrum form. If the 
rotation wavenumber Jt n lies in the inertial subrange, then for wavenumbers less 
than k n the turbulence motions are affected by rotation and the energy spectrum 
slope is modified. The present findings provide a new insight into the nature of the 
rotation effects on turbulence and, needless to say, their confirmation by (numerical) 
experiments would be desirable. It may, however, be difficult to experimentally 
distinguish the change in the spectral slope around the rotation wavenumber. The 
energy decay laws inferred in CSMC and the present results suggest a modification 
of the e model equation and eddy viscosity in k-e models. This is a subject of the 

following note. 
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2.2 A note on the eddy viscosity in rotating turbulence 

A suggested generalization of the expression for the eddy viscosity (in the consti- 
tutive relation u,uj = -2 vtS %} 4- \q 2 Sij) in rotating turbulence derived by Zeman 
(1994b) (Eq. (4) in preceding Section 2.1) is 


1 + — eijkRij)r 2 ' ' 

Here, vy 0 is the appropriate eddy viscosity for nonrotating flow; otherwise the 
notation is as in Section 2.1, i.e., q 2 = u]u] is twice turbulent kinetic energy (TKE), 

T ~ 2 / € ’ Rij = 5(^*1 i ~ Uj,i) the rotation tensor. The optimal value of the 
numerical constant c 5 was found to be c 5 « 0.1. 

In homogenous rotating turbulence with shear S = dU x /dx 2 = 2 R 12 and with 
the reference frame rotation Cl } = Cl 3 , (11) reduces to 


UT — 


VTo 

1 + c 5 Cl(Cl - 5/2)r 2 ' 


( 12 ) 


For a given value of the rotation-free viscosity vto and assumed constant value 
of the normalized shear Sr = Sq 2 /e , the eddy viscosity u-f is solely a function of 
the ratio Cl/ S and reaches maximum when Cl/S = 0.25 in agreement with linear 
stability analysis. The function vt{CI/S) is symmetric about Sl/S = 0.25 and falls 
off rapidly with increasing departure from 0.25. With the (tested) value of Sr = 12, 
ut decreases by a factor of 14 as Cl/S changes from 0.25 to 0.2510.75. 

Apart from the eddy viscosity, the rotation also affects the Kolmogorov energy 
cascade and therefore the rate of dissipation. The author proposed a correction to 
the e equation to represent this effect (reported also in Hadid, Mansour & Zeman 
1994). In the case of purely decaying turbulence subjected to rotation, the modified 
e equation is 


de 

dt 



(13) 


where /? is now a function of the rotation parameter u = |fl|r, in the following way 


/? = 3.7 + 


5 lo 2 
3 1 + uj 2 


(14) 


In this formulation, (13) satisfies the decay law dq 2 /dt oc t~ n so that the energy 
decay exponent n is 1.2 when w = 0 and n = 0.6 when u » 1. The latter value 
is based on the asymptotic decay of rotating turbulence inferred from the scaling 
analysis and LES results of Squires et al. (1993) (when the energy spectrum E(k) 
at the large scale end behaves as £ oc k*). The form of the function (14) has been 
based on the analysis of Zeman (1994b). The model-experiment comparison for 
rotating decaying turbulence using (13) and (14) is shown in Fig. 3. The data are 
from the experiment of Wigeland & Nagib (1978). 


L3 


FIGURE 3. Decay of rotating turbulence; model-experiment comparison. Data 

points are from Wigeland & Nagib (1978). i) ro = 0.12 ( ),0.47 ( ), 70.0 

( ); n To = 0.12 (■ ), 0.47 (* ), and 70.0 (• ). 



5, 

FIGURE 4. Shear-driven turbulence with rotation: k-e model prediction (with 
modified viscosity in (12)); cross-hatched areas represent roughly the DNS results 
of Bardina et a l. (1985). ft/S = 0.0 ( ), 0.25 ( ), 0.5 ( ), and -0.5 
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In the presence of shear 5 = dU\/dx 2 , the relevant k — e model equations are 


l*L = P 

2 dt 




(15) 


and (13) changes to 

^ = -/?(£ -0.75 P)^. (13') 

The TKE production rate P = 4i 'rSijSij = utS 2 is determined with the aid of 
the modified eddy viscosity in (12); the rotation parameter in (14) now has to 
include the contribution due to the presence of shear, i.e. u? 2 = (fi* — *ijkRij) 2 T 2 = 
(fi — S/2) 2 t 2 . The results of comparison between the model represented by (13’) 
and (15) and the DNS results of Bardina et al (1985) are presented in Fig. 4. The 
trend in the turbulence evolution with varying (l/S is apparently predicted although 
the k — e model is incapable of predicting the rapid distortion regime during the 
initial development when St < 1. It is noted, however, that the majority of the 
Reynolds stress closure models are incapable of reproducing the il effect on shear 
turbulence, particularly for the case of maximum amplification when f 1/S = 0.25. 
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